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Abstract. We investigate the mid-infrared radiation of galaxies that are powered by a starburst or by an AGN. 
For this end, we compare the spectra obtained at different spatial scales in a sample of infrared bright galaxies. 
ISO observations which include emission of the nucleus as well as most of the host galaxy are compared with 
TIMMI2 spectra of the nuclear region. We find that ISO spectra are generally dominated by strong PAH bands. 
However, this is no longer true when inspecting the mid-infrared emission of the pure nucleus. Here PAH emission 
is detected in starbursts whereas it is significantly reduced or completely absent in AGNs. A physical explanation 
of these new observational results is presented by examining the temperature fluctuation of a PAH after interaction 
with a photon. It turns out that the hardness of the radiation field is a key parameter for quantifying the photo- 
destruction of small grains. Our theoretical study predicts PAH evaporation in soft X-ray environments. Radiative 
transfer calculations of clumpy starbursts and AGN corroborate the observational fact that PAH emission is 
connected to starburst activity whereas PAHs are destroyed near an AGN. The radiative transfer models predict 
for starbursts a much larger mid-infrared size than for AGN. This is confirmed by our TIMMI2 acquisition images: 
We find that the mid-infrared emission of Seyferts is dominated by a compact core while most of the starbursts 
are spatially resolved. 
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1. Introduction 

Luminous infrared galaxies in the local universe with in- 
, frared luminosities exceeding 10^^ Lq were discovered in 
■ the IRAS all sky survey (Soifer et al. 1987). Responsible 
for their huge luminosity may be a combination of star- 
burst and AGN activity. In both cases, the nuclear emis- 
sion is fuelled by an enormous concentration of gas and 
dust in the central region (Sanders & Mirabel 1996). 

Luminous infrared galaxies are also found to be a ma- 
jor component in the early universe (z > 1). Deep mid- 
infrared surveys with the infrared space observatory, ISO, 
could identify individual objects which account for most 
of the cosmic infrared background light. It is argued that 
the emission of most of the distant infrared galaxies comes 
mainly from starbursts (Elbaz et al. 2002). This hypoth- 
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esis may also be tested by inspecting the spectra of lo- 
cal galaxies. ISO spectra of local galaxies (Rigopoulou 
et al. 1999, Siebenmorgen et al. 1999, Clavel et al. 2000, 
Laureijs et al. 2000, Laurent et al. 2000) usually show a 
set of infrared emission bands which are attributed to the 
emission of polycyclic aromatic hydrocarbons (PAH) . 

The question as to which source, on average, is respon- 
sible for the luminosity of the cosmological background 
is still a matter of intense debate and, at the moment, 
different answers are given for different wavelength re- 
gions. The X-ray background could be partly resolved 
by ROSAT into individual components and Hasinger et 
al. (1998) showed that at X-rays AGNs make the greatest 
contribution. At infrared wavelengths, on the other hand, 
starbursts seem to dominate, as infered from numerous 
ISO observations (Genzel & Cesarsky 2000). 

We want to re-address the question whether in- 
frared bright galaxies are predominantly powered by dust- 
enshrouded AGN or starbursts. As a new observational 
technique, we employ mid-infrared sub-arcsec imaging 
and mid-infrared spectroscopy at various spatial scales. 
(For previous ground based mid-infrared spectroscopy see 
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for example: Smith et al. 1989, Roche et al. 1991 and 
Soifer et al. 2002). The mid-infrared wavelength region 
is of particular interest as the optical depth is usually suf- 
ficiently low to penetrate the dust enshrouded nucleus. 
Modern mid-infrared cameras have adequate spatial reso- 
lution to estimate the size of the emitting regions which is 
an important parameter to constrain the nature of the un- 
derlying source: At comparable distances of the sources a 
resolved mid-infrared emission points towards a starburst 
galaxy while an unresolved nucleus indicates an AGN or 
an extreme compact star cluster. As we will demonstrate 
in model calculations, infrared emission bands are a diag- 
nostic tool to derive the activity type. The ratio of AGN 
versus starburst luminosity of local galaxies may be quan- 
tified by mid infrared spectra observed at different spatial 
resolutions. 

The paper is structured as follows: In Sect. 2 we in- 
troduce the sample of infrared bright galaxies observed 
with ISO and TIMMI2. The observation are summarised 
in Sect. 3, the data analysis procedure is explained in 
Sect. 4. Observational results of individual galaxies are 
given in Sect. 5. Sect. 6 and 7 presents radiative transfer 
calculations, which also illustrate the effect of the optical 
depth on the SED, and a discussion of the heating and 
evaporation of PAH in starburst and AGN environments. 

2. The sample 

The ISO archive contains more than 250 galaxies that 
have been spectroscopically observed in the mid-infrared. 
We have selected 22 of these galaxies for observation with 
TIMMI2 to compare the large aperture ISO spectra to 
the narrow slit TIMMI2 spectra. The selected galaxies all 
have 10 /im ISO flux densities above ~100mJy and are 
visible from La Silla (—29° South). Basic properties of the 
sample are presented in Table [3 The flux criterion was 
prompted by the sensitivity of TIMMI2 in spectroscopic 
mode which is ^lOOmJy/lOcr in 1 hour. The sample is 
heterogenous with respect to the galaxy type: it contains 
six starbursts (SB), one HII, two Syl and thirteen Sy2. 
The distances range from 3.5-176 Mpc and the luminosi- 
ties from 0.7 - 120 x IQI^Lq. Our TIMMI2 sample also 
contains sources lacking an ISO spectrum. These sources 
were selected from the Markarian catalogue (Markarian et 
al. 1989) to fill gaps in our primary observing schedule. 

3. Observations 

The N band imaging and spectroscopic observations were 
made with the TIMMI2 mid-infrared camera (Reimann et 
al. 2000, Kaufl et al. 2003) at the ESO 3.6m telescope. The 
camera uses a 240 x 320 pixel Raytheon Si:As array. In 
imaging we selected the 0.2" pixel scale providing a total 
field of view of 48" x 64". For long-slit spectroscopy the 
resolving power is AA/A ^ 400 and we used the 1.2" or 3" 
slit and 0.45" pixel field of view. Table |2 reports for each 
galaxy: observing date, total on-source integration time, 
filter and calibration standard. The filter is designated by 



Table 1. A few parameters of our IRG sample. Velocities 
(col. 3) and distances (col. 4) are from NED. Col. 5 gives 
the linear size of 1" at the distance D, and col. 6 the 
infrared luminosity. 



Name 


TvDe 


V 


D 


1" 


L 






km/s 


Mpc 


DC 




V^t^±± Let LIL Li.o a\. 




547 


3.5 


17 


0.68 




.Sv2 n 


439 


4 


19 


1.4 


IC4329a 


Sv 1 2 


4813 


65 


315 


7.6 


IRAS05189— 2524 


Sy 2.0 


12760 


176 


852 


120 


TR Asnsnn?— fifinn 


HII 


12324 


170 


822 


34 


M83 


SB 


516 


4 


19 


0.65 


Mrk.^OQ 

J.VX1 I\.\J\J \J 


Sv 1 2 


10312 


141 


684 


14 


Mrkl093 


SB 


4441 


60 


290 


9.3 


MrkMfifi 

J.VX1 iV X t: \J yj 


SB 


1327 


18 


87 


0.7 


NGC1068 


Sv2 

uy ^.Kj 


1148 


14 


68 


19 




Sv 1 8 


1636 


22 


106 


8.3 


NGC1386 


Sv2 

uy ^.Kj 


868 


18 


87 


0.71 


NGC2966 


SB 


2044 


27 


131 


1.1 


NGC3256 


SB 


2738 


37 


178 


35 


NGC4388 


Sy2.0 


2524 


34 


164 


1.2 


NGC5506 


Syl.9 


1853 


25 


120 


2.3 


NGC5643 


Sy2.0 


1199 


16 


78 


1.8 


NGC6000 


SB 


2193 


29 


141 


9.0 


NGC6240 


Sy2.0 


7339 


100 


483 


60 


NGC7552 


SB 


1585 


21 


102 


9.8 


NGC7582 


Sy2.0 


1575 


21 


102 


6.3 


NGC7674 


Sy2.0 


8671 


118 


573 


31 


PKS2048-57 


Sy2.0 


3402 


46 


222 


6.1 



its central wavelength and the entry grating refers to 8-13 
fim spectroscopy. 

Observations were performed in chopping and nod- 
ding mode with 12" amplitude each with is a compro- 
mise between sensitivity and the ability to measure ex- 
tended structures. To avoid saturation of the detector by 
the high ambient photon background the elementary inte- 
gration times were typically 20 ms for imaging and 40 ms 
for spectroscopy. Elementary images were coadded in real 
time during the chopping cycles. Chopping was in North- 
South direction at a frequency of ~ 3Hz. After an interval 
of typically less than a few minutes the telescope was nod- 
ded in order to cancel the telescope emission and residuals 
in the sky subtraction. After coaddition of all chopping 
and nodding pairs there are two negative and one posi- 
tive beam (slit) on the image which we combine applying 
a shift-and-add procedure. We note that TIMMI2 is not 
capable of measuring low surface brightness emission that 
is extended by more than half the chopping or nodding 
amplitude. Consequently, part of the diffuse emission is 
cancelled by the observing technique. 
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Target 


Date 


Time 


Filter 


Calibrator 




(2002) 


(sec) 


(Mm) 




Centaurus A 


Jan. 30 


645 


10.4 


HD81420 




Jan. 30 


3225 


grating 


HD81420 


Circinus 


Jan. 30 


268 


10.4 


HD81420 




Jan. 30 


967 


grating 


HD81420 


IC4329A 


Jan. 29 


403 


10.4 


HD81420 




Jan. 29 


2419 


grating 


HD81420 


IRAS05189-2524 


Feb. 01 


361 


8.6 


HD32887 




Feb. 01 


2257 


grating 


HD32887 




25.01.03 


600 


10.4 


HD123139 




25.01.03 


600 


11.9 


HD123139 


IRAS08007-6600 


Feb. 01 


361 


8.6 


HD32887 




Feb. 01 


469 


11.9 


HD32887 




Feb. 01 


1935 


grating 


HD32887 




25.01.03 


360 


11.9 


HD123139 


M83 


Jan. 30 


322 


10.4 


HD81420 




Jan. 30 


301 


11.9 


HD81420 




Jan. 30 


2419 


grating 


HD81420 




Jan. 30 


322 


10.4 


HD81420 


Mrk509 


Aug. 14 


619 


11.9 


HD196171 




Aug. 14 


619 


11.9 


HD196171 




Aug. 14 


2580 


grating 


HD196171 


Mrkl093 


Feb. 02 


602 


11.9 


HD32887 




Feb. 02 


2903 


grating 


HD32887 


Mrkl466 


Feb. 04 


903 


8.6 


HD123139 




Feb. 02 


903 


11.9 


HD123139 




Feb. 02 


2419 


grating 


HD123139 


NGC1068 


Jan. 29 


268 


10.4 


HD32887 




Jan. 29 


1209 


grating 


HD32887 




Jan. 31 


150 


8.6 


HD32887 




Jan. 31 


604 


grating 


HD32887 


NGC1365 


Jan. 30 


645 


10.4 


HD32887 




Jan. 30 


3225 


grating 


HD32887 




Jan. 31 


150 


8.6 


HD32887 




Jan. 31 


1935 


grating 


HD32887 


NGC1386 


Aug. 15 


619 


11.9 


HD1522 




Aug. 15 


2580 


grating 


HD1522 


NGC2966 


Feb. 03 


602 


11.9 


HD123139 


NGC3256 


Jan. 31 


180 


8.6 


HD123139 




Jan. 29 


537 


10.4 


HD81420 




Jan. 29 


3628 


grating 


HD81420 




Jan. 31 


3064 


grating 


HD123139 


NGC4388 


Feb. 03 


602 


11.9 


HD123139 




Feb. 03 


2419 


grating 


HD123139 


NGC5506 


Feb. 02 


602 


11.9 


HD123139 




Feb. 02 


1612 


grating 


HD123139 


NGC5643 


Feb. 02 


602 


11.9 


HD123139 




Feb. 02 


1612 


grating 


HD123139 




Aug. 14 


602 


11.9 


HD123139 




Aug. 14 


1612 


grating 


HD123139 


NGC6000 


Jan. 29 


215 


10.4 


HD123139 




Feb. 01 


301 


8.6 


HD123139 




Feb. 01 


2096 


grating 


HD123139 


NGC6240 


Aug. 13 


619 


11.9 


HD123139 




Aug. 15 


619 


11.9 


HD 133774 




Aug. 13 


3096 


grating 


HD123139 




Aug. 15 


2580 


grating 


HD133774 



Target 


Date 


Time 


Filter 


Calibrator 




(2002) 


(sec) 


(/im) 




NGC7552 


Aug. 15 


929 


11.9 


HD1522 




Aug. 15 


2580 


grating 


HD169916 




Aug. 15 


3096 


grating 


HD1522 


NGC7582 


Aug. 14 


619 


11.9 


HD178345 




Aug. 14 


2580 


grating 


HD178345 




Aug. 14 


2709 


grating 


HD1522 


NGC7674 


Aug. 15 


619 


11.9 


HD 178345 


PKS2048-57 


Aug. 14 


619 


11.9 


HD178345 




Aug. 14 


2460 


grating 


HD178345 



TIMMI2 follow-up observations of ISO galaxies were 
obtained during January 28-February 3, August 13-14, 
2002 and 4 photometric imaging observations on January 
25, 2003. The observations of the targets were interleaved 
with those of nearby mid-infrared calibration standard 
stars, typically 4-5 of which were measured each night. 
Model spectra for the standards have an absolute photo- 
metric uncertainty of a few percent (Cohen et al. 1999). 
However, for stars of spectral type later than F there 
might be fundamental molecular lines present in the pho- 
tospheric spectra, such as Si-0 in the 9/xm region (Heras 
et al. 1998). These bands are not considered by the present 
photospheric models. Consequently, one must be care- 
ful not to over-interpret faint structures at these wave- 
lengths. Furthermore close to the blue and red cut-off fre- 
quencies of the grating the terrestrial background varies 
strongly with ambient conditions and atmospheric correc- 
tions are more difficult. The same holds for wavelength 
regions of strong terrestrial lines near 9.58/im, 11.73/um 
and 12.55/xm. 

The N band imaging of the stars was performed in 

the same filter as the target and used for photometric 
flux conversion from photon count rates into astronomical 
units (mJy) and to establish the point spread function 
(PSF). The spectra of the stars were used for atmospheric 
corrections and for flux calibration of the target spectrum. 

4. Data analysis 

4.1. Mid-infrared imaging 

We apply multi-aperture photometry on the shift-and- 
add processed image, which we refer to as raw or original 
image. For each aperture centred on the brightest pixel of 
the source we estimate the background as the mean flux 
derived in a 6 pixel wide annulus which is put 6 pixels away 
from the outer source aperture. The thus estimated aper- 
ture fluxes, F{r), were plotted against the aperture radii, 
r, to provide the growth curve. The flux where the growth 
curve of the aperture photometry flattens was identifled 
interactively and is taken as the total flux density. The 
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Table 3. Mid-infrared flux densities of TIMMI2 observa- 
tions. Photometric uncertainty ~10%. 



N 8-1116 


Ff8 61 


FflO 4") 


Ffll 91 






fJv) 


fJv) 


C6iitcLiiriis A 




0.65 








4.10 




Tr!4S2QA 




0.64 




TR ASn.'il SQ— 2524 


0.42 


0.42 


0.57 


TR Asnsnn?— fifinn 


0.21 




0.23 


Mrk^OQ 






0.22 


MrklDQS 

i.V±l JV W t-/ o 






0.16 


Mrkl4fifi 


0.09 




0.11 


M83 




0.11 


0.23 


NGC1068 


14.8 


17.6 






0.40 


0.46 










0.40 








<0.03 




0.45 


0.39 




NOC4S88 

J. i v_> ^ (J i_i i_) 






0.29 


NGC5506 


- 


- 


1.06 


NGC5643 






0.31 


NGC6000 


0.21 


0.20 




NGC6240 






0.28 


NGC7552 






1.87 


NGC7582 






0.69 


NGC7674 






0.26 


PKS2048-57 






0.92 



absolute photometric uncertainty is ^ 10%. This error 
estimate is based on internal consistency of the calibra- 
tion observations and monitoring of photometric calibra- 
tion standard stars. We note that part of this error is due 
to large scale (> 10") sensitivity gradients on the detec- 
tor array. Up to now no suitable flat flelding procedure 
has been established to further reduce the absolute pho- 
tometric error. 

The growth curve is used to extract spatial informa- 
tion. For example, if one calculates the growth curve of an 
Airy function up to its first dark ring then the aperture 
radius, r^Q%^ where the flux is at 50% of its maximum 
value is half of the full width at half maximum (FWHM). 
Still the diffraction pattern of the telescope and the seeing 
has to be subtracted. If for an extended target one derives 
a 50% flux radius rt^Q% and for a PSF a radius r^^y^, a size 
estimate of the source extension may be approximated by 

(''50%)^ = (''soro)^ - ('"50%)^- The size estimate by the 
growth curve depends on the total flux or a normalisation 
at large radii. Another diagnostic of the source extension 
is given by the differential flux profile, which is a func- 
tion F{ri+i) — F{ri), and as well difiicult to normalise, 
this time at small radii. We take the sum of the aperture 
fluxes at r < r^Q%- With such normalisation the differen- 
tial flux distribution of an extended target should be less 
steep than for the PSF. 



The spatial resolution of the images are at best limited 
by the diffraction pattern of the 3.6m telescope, which 
provides a FWHM of 0.7" at lO/xm. In oder to increase the 
spatial resolution we applied a deconvolution procedure. 
We use the maximum entropy multi-resolution wavelet 
decomposition by Starck & Murtagh (2002). The gradient 
algorithm has been chosen with 4 wavelet scales assuming 
Gaussian distribution of the noise and we applied a 5a 
threshold, where a is the noise as estimated on the object 
free part of the image. The method provides a residual 
map, R, which is a convolution of the deconvolved image, 
D, with the PSF, P, and subtracted from the original 
image, O; symbolically: R = D * P — O. The iteration of 
the deconvolution is stopped so that the noise properties 
of the residual map is at similar level as the noise of the 
original image. However, for bright point-like sources with 
little extension we note that the residual map is polluted 
by the very bright central peak. Residuals are typically less 
than 5% of the peak flux of the deconvolved image. They 
are caused by finite sampling of a 5-type function. From 
the growth curve of the deconvolved image we calculate 
its 50% flux radius, r^Q%- The residual map is used for 
error propagation 

To study the stability of the PSF and its influence on 
the deconvolution procedure, we calculate for each PSF 
the deconvolution and growth curves with any other PSF 
available of the same fllter. For the deconvolved PSF im- 
ages we flnd a mean size of {r^^yj = 0.098" ± 0.046". 
Size estimates from the growth curves as measured on 
the raw images give a mean of {r?^Q%) = 0.224" ± 0.074", 
which already closely matches the static aberration of the 
telescope. Altogether we verified that seeing in the mid 
IR gives only a small contribution to the PSF at the 
3.6m. Nevertheless for deconvolution we apply the mea- 
sured PSF matched closely in time to the target obser- 
vation. A source is regarded to be unresolved if its ra- 
dius measured on the raw and deconvolved image is be- 
low (^50%) + 3cr(r5o%). Therefore, resolved sources have 
^^50% > O-^^" ^50% > 0-44", respectively. For galactic 
nuclei fulfilling these criteria, the differential source profile 
also shows signatures of extended emission. 

Observed mid-infrared flux densities for the galactic 
nuclei are given in Table 3. For 10 out of 22 galaxies the nu- 
cleus is resolved. For these we show in Fig.lTHTTlraw. PSF, 
deconvolved and residual maps and present growth curves 
as well as differential flux distributions. Absolute levels 
of the contours, shown on top of the gray scale images, 
can be derived together with Table 4, where we specify 
the peak flux densities of the raw and deconvolved images 
and the residual map. For the residual map we specify the 
la RMS noise level. The la RMS noise of the raw im- 
age is given in col. 3 of Table 5. It is an indicator of the 
sensitivity limit which we reached with TIMMI2. 

Table 5 lists the central wavelength of the TIMMI2 
filter (col. 2) together with the 50% flux radius, r^Q%, as 
derived from the raw (col. 4) and deconvolved (col. 5) 
image. The latter radius is our estimate of the source size 
(col. 6). 
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Table 4. Peak surface brightness and RMS of nuclei with 
gsctcnded mid infrared emission. 
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Fig. 1. NGC3256 at 10.4/im. On the raw (top-left), PSF 
(top-middle) and deconvolved image (bottom-left) con- 
tours are in black at 10%, 30% and 50% and in white at 
75%, 90% and 95% levels of the peak flux; dotted contours 
are at —la and —3a and dashed contours at +la and +3a, 
respectively. On the residual map (bottom-middle) black 
contours are at —98%, —90%, —50%, white contours at 
50%, 90%, 98% of the peak. The peak fluxes and a values 
are specified in Table 4. For all gray scale images North is 
up and East is left. Growth curves (top-right) and differ- 
ential flux distributions (bottom-right) are shown for raw 
(dashed), PSF (dotted) and deconvolved image (full line). 
Error bars are la RMS. 

Fig. 2. NGC7552 at 11.9 /im, see Fig.[T] On the raw image 
(top-left), the PSF and the deconvolved image contours 
are set in black at 8%, 20% and 30% and in white at 50%, 
75%, and 95% levels of the peak flux 

Fig. 3. Mrkl093 at 11.9 /xm, see Fig.H 

Fig. 4. Circinus at 10.4 /im, see Fig. ^ 

Fig. 5. Centaurus A at 10.4/im, see Fig.^ 

Fig. 6. NGC1068 at 10.4/im, see Fig.[T] 

Fig. 7. NGC1386 at 11.9/im, see Fig.[T] 

Fig. 8. NGC7582 at 11.9 /im, see Fig.[T] On the raw image, 
the PSF and the deconvolved image contours are set in 
black at 14%, 20% and 50% and in white at 75%, 90%, 
and 95% levels of the peak flux, respectively. 

Fig. 9. Mrk509 at 11.9 /im, see Fig. [2 
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Name 


Raw 


Deconvolved 




Residual 






Peak 


Peak 


Peak 


RMS 






mjy 


mJy 


mJy 


mJy 




1 


arcscc 


arcscc 


arcsec 


5 

arcsec 




Centaurus A 


881 


2301 


100 


17 




Circinus 


3902 


14806 


167 


31 




Mrk509 


637 


10 


2.6 


0.4 




Mrkl093 


52 


10 


7.5 


1.8 




NGC1068 


10078 


31196 


371 


80 




NGC1386 


512 


136 


33 


4.7 




NGC3256 


296 


211 


35 


7.3 




NGC6240 


280 


245 


56 


11 


NGC7552 


997 


277 


110 


10 


NGC7582 


275 


91 


12 


2.3 


NGC7674 


287 


123 


16 


4.0 



Table 5. Mid-infrared sizes of galactic nuclei. 



1 


2 


3 


4 


5 


6 


Name 


A 


a" 






'^50% 






mJy 

arcscc^ 


ff 


II 


pc 


Centaurus A 


10.4 


14 


<0.37 


<0.20 


<3 


Circinus 


10.4 


22 


0.46 


<0.20 


<4 


IC4329a 


10.4 


19 


<0.28 


<0.21 


<66 


IRAS05189-2524 


8.6 


18 


<0.24 


<0.20 


<170 


IRAS08007-6600 


8.6 


66 


<0.23 


<0.27 


<222 


Mrk509 


11.9 


10 


1.2 


2.4 


1642 


Mrkl093 


11.9 


13 


1.9 


1.69 


490 


NGC1068 


8.6 


48 


0.53 


0.27 


18 


NGC1068 


10.4 


18 


0.58 


0.35 


24 


NGC1365 


8.6 


64 


<0.28 


<0.35 


<37 


NGC1386 


11.9 


42 


0.50 


0.53 


46 


NGC3256 


10.4 


28 


1.14 


0.88 


157 


NGC4388 


11.9 


33 


<0.30 


<0.36 


< 59 


NGC5506 


11.9 


22 


<0.30 


<0.20 


<24 


NGC6240 


11.9 


31 


1.11 


0.47 


227 


NGC7552 


11.9 


82 


2.79 


2.65 


270 


NGC7582 


11.9 


17 


1.18 


0.87 


89 


NGC7674 


11.9 


22 


0.6 


0.46 


264 


PKS2048-57 


11.9 


15 


<0.2 


<0.23 


<51 



la RMS noise measured on the object-free part of the raw 
image. 

° Mid-infrared radius where the flux is at 50% of the total as 
derived from the growth curves of the object (raw) and PSF 
images. 

Mid-infrared radius where the flux is at 50% of the total as 
derived from the growth curve of the deconvolved image. 
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Fig. 10. NGC6240 at 11.9 ^m, see Fig.ffl 
Fig. 11. NGC7674 at 11.9/xm, see Fig.[T] 

4.2. TIM M 12 spectroscopy 

The coaddition of all chopping and nodding pairs of raw 
frames in TIMMI2 spectroscopic observing mode gives two 
negative and one positive long-slit spectrum on the im- 
age. We developed an optimal extraction procedure for 
TIMMI2 spectra similar to Horne (1986): A source profile 
is extracted by collapsing the image along the dispersion 
direction and applying a median filter. Extension of the 
sources are considered to a 2a cut-off, where a is calcu- 
lated for each column in cross-dispersion direction and on 
the source free part of the array. Preliminary wavelength 
calibration is performed using tables provided by the in- 
strument team. This calibration is fine tuned to sub-pixel 
accuracy by applying a linear wavelength shift. The shift 
is measured as an offset to atmospheric lines. The proce- 
dure is applied to the target and the calibration standard 
star. A division of both spectra removes the telluric lines. 
To minimise residuals of the sky line cancellation, we ob- 
served standard and target at similar airmasses. From the 
standard star spectrum we calculate the conversion factor 
from technical (ADU/s) into astronomical (Jy) units. The 
calibrated spectra are finally high-frequency noise filtered 
by the method described by Starck et al. (1997). Statistical 
errors are computed from unfiltered data. Because not all 
detector channels were always fully operational, the spec- 
tra of M83, CenA, Circinus, IC4329A, and NGC1068 could 
only be extracted long-ward of 9/xm. 

4.3. ISO spectroscopy 

Whenever possible, we compare mid-infrared spectra at 
different spatial scales: high-resolution long-slit TIMMI2 
spectroscopy (1" — 3") and ISO spectra of much lower 
resolution (24" x 24" with ISOPHT, Lemke et al. 1996; 
14" X 20" with ISOSWS in the N band, de Graauw et 
al. 1996; and ISOCAM spectra which refer, unless other- 
wise stated, to the total emission of the galaxy, Cesarsky et 
al. 1996). The diffraction limits of the telescopes at 10/zm 
are: ~ 0.7" at the ESQ 3.6m and ~ 8" with ISO. 

For the brightest objects of our sample ISOSWS 
grating or line scans are reported by Thornley et al. 
(2000), Sturm et al. (2002) and Verma et al. (2003); 
ISOCAM spectra by Laurent et al. (2000), and ISOPHT 
spectra by Clavel et al. (2000), Laureijs et al. (2000), 
Rigopoulou et al. (1999), Siebenmorgen et al. (1999) 
and Spoon et al. (2002). For the remaining galaxies 
(CenA, Circinus, IC4329A, IRAS08007-6600, Mrk509, 
NGC1068, NGC1365, NGC1386, NGC3256, NGC4388, 
NGC5506, NGC5643, NGC6000, NGC6240, NGC7552, 
NGC7582) we extracted ISOPHT raw data from the ISO 
archive and reduced the spectra with standard routines 
of the ISOPHT interactive analysis system (PIA). Most 



Table 6. Important features and lines between 8 to 13/im. 



Name 


A" 


I' 


Mode 




^m 


eV 




PAH 


8.6 




C-H 


Ar HI 


8.991 


27.6 




Mb VII 


9.009 


186.5 




Fe VII 


9.537 


99.1 


- 


Ha 0-0 S(3) 


9.660 






Silicate 


9.7 




Si-0 


S IV 


10.511 


34.8 




PAH 


11.3 




C-H 


Ne II 


12.814 


21.6 




PAH 


12.7 




C-H 



" Rest wavelength. 

^ lonisation energy with respect to the next lower ionisation 
level. 

of the ISOPHT spectra appear for the first time in the 
refereed literature. 

5. Results 

Here we present our high resolution mid-infrared imag- 
ing and spectroscopy data using TIMMI2 together with 
the low resolution ISO spectra. The data of the individual 
objects are briefly discussed and whenever necessary com- 
pared to observations available at other wavelengths. We 
will first describe the starbursts and then the Seyferts. 

A variety of broad band dust features and fine struc- 
ture lines fall within the N band; they are summarised in 
Table 6. Gas emission lines imply either the presence of 
ionising shocks (Contini & Contini 2003) or a hard pho- 
ton flux. For ionisation potential below ^50cV, photo- 
ionisation by hot stars is the most likely explanation for 
the excitation mechanism. In the 8-13/im spectral range 
there are no lines which require very hard radiation only 
produced by AGNs (> 200 eV). At the resolving power 
of TIMMI2 (A A/ A 400) the hues are generally unre- 
solved. The [Ar III] and [Mg VII] line occur at the same 
wavelength but their excitation requirements are very dif- 
ferent. The [Fe VII] line may be blended by the Ozone 
band, [S IV] is in a clean atmospheric region. The PAH 
band emission at 12.7/im is close to the [Ne II] at 12.8/im. 

However, as we will demonstrate by model calculations 
in Sect.[7| detection or non-detection of PAH features can 
serve as a diagnostic tool to distinguish the dominating ex- 
citation within the nucleus: starbursts or AGN; although 
the dividing line is not sharp. The mid IR spectra of spi- 
rals show dust continuum and strong PAH bands. Close to 
massive stars there is additional hot dust continuum emis- 
sion and some photo-destruction of PAH (Siebenmorgen 
1993). So that starburst galaxies still show PAH bands but 
at lower PAH band-to-continuum ratios than observed in 
spirals. In the harsh environment close to an AGN one ex- 
pects much lower PAH band-to-continuum ratios than in 
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starbursts and most often PAH bands are absent. In the 
immediate environment of the hard X-ray and intense UV 
radiation of an AGN the PAH get destroyed. In addition 
significant emission by hot dust from the inner torus wall 
of the unified model is expected which further decreases 
the PAH band-to-continuum ratios. 

The silicate band if present in the galactic nuclei is 
seen so far only in absorption. From its depth one can 
estimate the visual extinction. We apply Ty ^ 18 ■ rg y^m. 
However, Tg^ifim is difficult to estimate without detailed 
model calculations of the PAH contribution to the mid- 
infrared continuum (Siebenmorgen et al. 2001). Whenever 
possible we interpolate a line to the spectra from ~ 8 to 
11 /im and so derive a crude estimate of rg.y^m and the 
total visual extinction. This procedure maybe uncertain 
within a factor 2. 

5.1. Starburst galaxies: imaging and spectroscopy 

5.1.1. NGC6000 

ISO observations and radiative transfer models for the 
dust in NGC 6000 are presented by Siebenmorgen et al. 
(1999). At 12/xm, the flux from IRAS is stronger than 
from ISOPHT which suggests that the mid IR source is 
extended beyond 24". TIMMI2 images at 8.6 and 10.4/^m 
are fuzzy, without a point source. Both spectra in Fig. 1121 
are dominated by PAH emission. In the TIMMI2 spec- 
trum, the 11.3/xm PAH band-to-continuum ratio is about 
two times smaller (Table 8). The [Ne II] line is de- 
tected, but the 12.7/im PAH band is seen only marginally. 
We estimate from the depth of the silicate absorption 
in the ISOPHT (TIMMI2) spectrum a visual extinction 
of ~ 25 mag (20 mag); both values are consistent with 
Ay — 29 mag inferred from the radiative transfer models. 

5.1.2. Mrkl466 

The signal-to-noise in the TIMMI2 image is low, hence 
we cannot determine the source size. The only feature in 
the TIMMI2 spectrum detected at high confidence is the 
ll.S^m PAH band (Fig.UHand Table 8). ISO spectra are 
not available for comparison. 

5.1.3. NGC3256 

The brightness in the TIMMI2 image increases towards 
the center. The emission is resolved and there is no point 
source in the core (Fig. P). Both ISOPHT and TIMMI2 
spectra are dominated by PAH emission (Fig. 112(1 . The 
11.3/^m PAH band-to-contimmm ratio is lower in the 
smaller aperture spectrum by about a factor 2 (Table 8). 
There is strong [Ne II] emission (Table 7). Due to the 
8.6/im PAH band and because the atmospheric ozone band 
at ~ 9.5/im has not been fully eliminated, it is difficult to 
say whether there is silicate absorption at all. 



Fig. 13. Mid-infrared flux densities (in Jy) of M83: 
Data presented are: ISOPHT (24" aperture, histogram), 
TIMMI2 spectroscopy (3", full line) with 3a errors 
(shadow) and TIMMI2 photometry with la error bars 
(Table 3). For comparison we show the spectrum of the 
Galactic Center (dashed). 

5.1.4. NGC7552 

The deconvolved 11.9^m TIMMI2 image at 0.45" angu- 
lar resolution (Fig. ^ shows four hot spots, equidistant 
from the center, that form a ringlike circumnuclear star- 
burst. Two of the hot spots have already been detected 
by Schinnerer et al. (1997) in their 10.5/im map with 0.8" 
resolution. These authors also present near infrared im- 
ages and spectral synthesis models consistent with the 
idea of a starburst ring. The slit during the TIMMI2 
grating observations was positioned at the brightest knot, 
^ 3" North of the centre; the [Ne II] line is visible as 
weU as the 11.3Aim PAH band (Table 7). The ISOPHT 
spectrum shows the [Ar II] line and strong PAH emis- 
sion. The 11.3/xm PAH band-to-continuum ratio is a fac- 
tor ~3 lower in the smaller aperture spectrum (Table 8). 
Silicate absorption is weak in both spectra, it corresponds 
to Av ^ 3 mag (Fis. T^ . 

5.1.5. Mrkl093 

The source at 11.9/im is extended in the East-West di- 
rection (Fig. 131 . TIMMI2 spectroscopy (Fig. [T^l) reveals 
the [Ne II] line and PAH bands at 8.6 and 11.3 ^m. ISO 
spectra are not available for comparison. 

5.1.6. IRAS08007-6600 

This is one of the most distant galaxy in our sample 
(170 Mpc). The TIMMI2 images at 8.6 ^m and 11.9 fim 
are unresolved. The energy distribution peaks near 60/im 
(Vader et al. 1993). The TIMMI2 spectrum does not show 
evidence for the presence of 8.6 /im PAH emission. This 
feature is, however, present in the larger aperture ISOPHT 
spectrum (Laureijs et al. 2000). The absence may be due 
to low signal-to-noise of our TIMMI2 spectrum. The sil- 
icate absorption is weak and indicates Ay ^ 6 mag from 
the central 200pc (Table 5). 

5.1.7. M83 

M83 is the nearest starburst galaxy in our sample. 
TIMMI2 images at 10.4 and 11.9/im show a fuzzy struc- 
ture without striking point sources. As typical for star- 
burst galaxies, the large aperture ISOPHT spectrum 
is dominated by PAH emission bands. In contrast, the 
TIMMI2 spectrum of the central 60 pc does not show any 
sign of PAH emission features and is dominated instead 
by a smooth continuum with a weak [Ne II] emission line 
(Fig. • The shape of the continuum is consistent with 
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Fig. 12. Mid-infrared flux densities (in Jy) of starburst galaxies: Data presented are: ISOPHT (24" aperture, his- 
togram), TIMMI2 spectroscopy (3", full line) with 3cr errors (shadow) and TIMMI2 photometry with la error bars 
(Table 3). 

Table 7. Integrated atomic hne fluxes and 3cr upper limits (10"^°W/cm2) of TIMMI2 and ISOSWS spectra. 



Name 


[S IV] 


[S IV] 


[Ne II] 


[Ne II] 


Aperture 


3" 


14" X 20" 


3" 


14" X 27" 


Centaurus A 


<2.4 


1.4'' 


22.6 ± 1.4 


22.1'' 


Circinus 


<3.6 


12.7'' 


6.6 ± 0.9 


90'' 


IRAS05189+2524 


<1.7 








IRAS08007-6600 


<1.2 








IC4329a 


1.8 ± 0.5 




1.6 ± 0.3 




M83 


<0.3 




1.7 ± 0.2 




Mrkl093 


<2.2 




3.0 ± 0.5 




Mrkl466 


<1.8 




<2.3 




Mrk509 


<1.3 




<1.3 


1.3'' 


NGC 1068 


30.1 ± 6.2 


58'' 


<14.2 


70'' 


NGC 1365 


<1.7 


2.6'' 


<1.8 


40.9'' 


NGC 1386 


<2.8 




<2.8 




NGC 3256 


<1.6 


0.9" 


28.5 ± 1.0 


89.2" 


NGC 4388 


3.4 ± 0.9 




5.4 ± 0.9 




NGC 5506 


3.5 ± 0.9 


5.4'' 


2.9 ± 0.9 


5.9'' 


NGC 5643 


<1.3 




1.6 ± 0.4 




NGC 6000 


<2.0 




15.1 ± 1.1 




NGC 6240 


<1.4 




5.0 ± 1.6 


17" 


NGC 7552 


<8.3 


0.3" 


49.0 ± 5.6 


68" 


NGC 7582 (1.2") 


<3.7 


1.8'' 


<8.8 


14.8'' 


NGC 7582 (3.0") 


<1.6 


1.8'' 


8.9 ± 1.2 


14.8'' 


PKS2048-57 


<2.7 




<3.6 


2.1'' 



" Thornley et al. (2000), '' Sturm et al (2002), " Verma et al. (2003) 



the profile of the silicate absorption feature as found to- 
wards the ISOSWS position Sgr A* in the Galactic center 
(Lutz 1999). Given their similarity, we speculate that the 
conditions in the nucleus of M 83 may be similar to those 
in the Galactic center. In the latter, PAHs are thought to 
be destroyed by the local intense radiation field (~ 10^-^ x 
ISRF) contributed by massive stars found in the Galactic 
center (Lutz 1999). Indeed far infrared spectroscopic im- 
ages of fine structures lines in M83 indicate a strong cen- 
tral starburst activity headed by 09 stars (Stacey et al. 
1998). So that this is our best example where PAH de- 
struction occurs by the strength of the radiation field 
rather than by its hardness. The latter process is discussed 
in Sect. 6. From the depth of the silicate feature we esti- 
mate a visual extinction '-~^20mag. 



5.2. Seyfert galaxies: imaging and spectroscopy 
5.2.1. Circinus 

Circinus is the closest galaxy of our sample that con- 
tains an obscured Seyfert nucleus. In the TIMMI2 10.4/im 
image (Fig. ^J, half the flux comes from an unresolved 
core, the other half from a l"-sized blob elongated North- 
South. At a comparable scale, resolved emission at 8.5/im 
is reported by Siebenmorgen et al. (1997). The photom- 
etry of ISOPHT (this work) and ISOSWS (Moorwood et 
al. 1996) give similar fluxes for A > 9^m, but ISOSWS is 
systematically weaker 20%-30% between 6 — 9/im. From 
the silicate absorption feature we estimate a visual extinc- 
tion Ay > 15 mag. 

In the large aperture ISOSWS spectrum, one sees the 
silicate absorption feature, strong line emission ([S IV], 
[Ne II] ) and PAH bands. In the TIMMI2 spectrum, [S IV] 
and [No II] arc not detected so they cannot be of nuclear 
origin. The TIMMI2 absolute fluxes agree well with Roche 
et al. (1991) made with a 5" aperture. 
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Fig. 14. Mid-infrared flux densities (in Jy) of Seyfert galaxies: Data presented are: ISOPHT (24" aperture, histogram), 
ISOCAM (total field or nuclear spectrum, dashed), ISOSWS (14- 20", dotted), UCL (5", diamonds), TIMMI2 spec- 
troscopy (3", full line) with 3cr errors (shadow) and TIMMI2 photometry with la error bars (Table 3). 



Moorwood & Oliver (1994) find that the central 6pc 
region is dominated by an AGN whereas star formation, 
as traced by Ha and Br7 recombination lines, occurs in 
Circinus in a (partial) ring of 200pc radius. This ring lies 
outside our map of Fig. 3. Model calculations of the dust 
re-processing of the AGN are presented by Siebenmorgen 
et al. (1997) and Ruiz et al. (2001). 

Generally speaking, when one observes a Syl galaxy 
with a large aperture, one has a free view to the inner wall 
of the torus and detects strong mid-infrared emission from 
hot (large) grains. The PAH band-to-continuum ratio is 
then relatively low. In a Sy2, like Circinus, the inner torus 
is obscured and hence, the hot grains are invisible and the 
PAH band-to-continuum ratio is higher. This picture is 
supported in the case of Circinus, but also for Syl and 
other Sy2 galaxies (see below). 

5.2.2. Centaurus A 

This is the closest example of an AGN in our sample. In 
the TIMMI2 image it shows an unresolved core of less than 
0.5" (Fig.jSJ. The raw image and the growth curve (Fig.lsj 
indicate weak extended emission over — 2" with a surface 
brightness ten times lower than at the peak. The ISOCAM 
(Laurent et al. 2000) and ISOPHT (this work) spectra are 
PAH dominated whereas the TIMMI2 spectrum is fea- 
tureless except for the presence of a strong [Ne II] line 
(Fig. I14|l . The visual extinction estimated from the sil- 
icate absorption detected in the ISOPHT and TIMMI2 
spectra is around 14 mag. 

5.2.3. NGC1068 

TIMMI2 resolves the nucleus of this prototype Sy2 galaxy 
(Fig. and our image of the core agrees with the long 
integrations by AUoin et al. (2000) using CAMIRAS and 
Bock et al. (2000) with MERLIN at Keck. However, we are 
not sensitive enough to confirm the individual hot cores 
found in those maps. PAHs are detected in the integrated 
ISOCAM (Laurent et al. 2000) and in our ISOPHT spec- 
trum, but in a smaller aperture they disappear (Fig. 
see also the ISOSWS spectrum of Lutz et al. 2000). The 
continuum is of similar strength for ISOSWS and TIMMI2 
indicating that the AGN is indeed the major source in the 
central Ikpc region. 

The ISOCAM spectrum of the central 9" by Laurent 
et al. (2000) is, like ISOSWS and TIMMI2, free of PAHs. 
The [Ne II] line is absent in TIMMI2, but strong in the 
ISOSWS spectrum, [S IV] is seen in both. The strength of 
the silicate absorption band increases with spatial resolu- 
tion; for the central 20pc, we estimate Ay > 5 mag. 



5.2.4. NGC5643 

This Sy2 galaxy shows high-excitation optical lines 
(Kinney et al. 1993). The emission in [O HI] and Ha is 
extended up to 20" (Schmitt et al. 1994). The 11.9^m 
TIMMI2 image is fuzzy, without a point source. The 
[Ne II] line is detected by TIMMI2 and for the [S IV] line 
we get an upper limit (Table 7). PAHs are strong in the 
ISOPHT aperture, whether they are seen with TIMMI2 
is unclear fFig. I14|l . 

5.2.5. NGC1386 

The TIMMI2 image of this Sy2 galaxy (Ferruit et al. 2000) 
has an unresolved core (Table 5) and an extended compo- 
nent (Fig. d). The ISOPHT spectrum shows PAH bands 
and tentative, weak silicate absorption. In the TIMMI2 
spectrum, (Fig. [Tl|l one sees the silicate feature, corre- 
sponding to Ay sa 30 mag, but no PAHs. 

5.2.6. NGC1365 

A face-on Syl galaxy with an unresolved core in the 
TIMMI2 image. Its mid-infrared emission is dominated by 
PAHs everywhere except in the very center. The band-to- 
continuum ratio decreases with decreasing aperture size 
(Fig.ini). PAH bands are not present in the TIMMI2 spec- 
trum (Table 8). [S IV] and [Ne II] emission is detected only 
with ISOSWS, but not with TIMMI2. The TIMMI2 con- 
tinuum is flat without evidence of silicate absorption. 

5.2.7. NGC4388 

A galaxy with a Sy2 nucleus where Ho & Ulvestad (2001) 
see at radio wavelengths two components, 1.9" (150 pc) 
apart, none of them coincides with the position of the op- 
tical peak. A hard X-ray peak is found at the position of 
the active nucleus (Iwasawa et al. 2003). As in the previ- 
ous examples, the TIMMI2 image has a single unresolved 
core (Table 5), and there are strong PAH bands in the 
ISOPHT, but none in the TIMMI2 spectrum (Fig. UH). 
The [S IV] and [Ne II] lines are detected with TIMMI2 
(Table 7). 

5.2.8. IRAS05189-2524 

In this most distant object of our sample, the central 
source in the TIMMI2 image is unresolved (Table 5). The 
same conclusion was reached by Soifer et al. (2000) us- 
ing MERLIN at Keck. ISOPHT sees weak PAH emission 
(Laureijs et al. 2000) and an absorption feature at Gfim 
attributed to water ice (Spoon et al. 2002). PAH bands 
(Fig. I14|l as well as gas emission lines are absent in the 
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Fig. 15. Mid-infrared flux densities (in Jy) of Seyfert galaxies: Data presented are: ISOPHT (24" aperture, histogram), 
TIMMI2 spectroscopy (3", full line and 1.2" dotted) with 3cr errors (shadow) and TIMMI2 photometry with la error 
bars (Table 3). 



TIMMI2 spectrum. From the silicate absorption, we esti- 
mate A\r ^ 12 mag. 

5.2.9. NGC7582 

The spectrum is a superposition of a Sy2 nucleus and a 
starburst (Storchi-Bergmann et al. 2000). The AGN is re- 
sponsible for X-ray emission over a few kpc (Levenson et 
al. 2000). Optical and near IR images (Regan & Mulchaey 
1999) indicate a complex dust morphology. There is a lane 
that runs from the southeast to an absorbing ring north of 
the chain of blue clumps; the TIMMI2 image reflects this 
structure (Fig.|Hl). PAH bands are detected with ISOPHT, 
but not [S IV] . The excellent seeing conditions during the 
TIMMI2 observations allowed us to use besides the 3" also 
the 1.2" slit. The spectra look quite similar at both resolu- 
tions fFig. [T5|l : [Ne II] is present, the PAH band is detected 
at 11.3/im but not at 8.6 /zm. From the silicate absorption 
we estimate Ay ^ 20 mag. The 11.3/Ltm PAH band-to- 
continuum ratio is in TIMMI2 lower than in ISOPHT by 
about a factor 6 (Table 8). 

5.2.10. NGC5506 

The optical spectrum of this Sy2 galaxy is comparable 
to what one finds in young starburst galaxies, however, 
the high excitation hues ([Ne V], [Fe VII], [He II]) man- 
date the presence of a non-thermal continuum (Oliva 
et al. 1999). The X-ray spectrum is highly absorbed 
(Pfefferkorn et al. 2001). In the 8-13/im wavelength range, 
[Ar III]/[Mg VII], [S IV] and [Ne II] emission lines are 
detected with ISOSWS, [S IV] and [Ne II] hues at simi- 
lar strengths with TIMMI2 (Table 7). Evidence of PAH 
emission is found in the ISOPHT spectrum at 6.2/xm 
(Rigopoulou et al. 1999), whereas TIMMI2 sees a fea- 
tureless spectrum. TIMMI2 and ISOPHT agree very well 
photometrically fFig. 1151) . We estimate an optical depth 
of over 15 mag towards the core. 

5.2.11. PKS2048-57 

In the near-infrared, the nuclear source of this object is 
very red (K-L = 2.1 mag) and heavily polarised (17%) 
which is explained by a non-thermal contribution (Hough 
et al. 1987). Its spectral energy distributions peaks at 
60/xm (Hessler & Vader 1995) and the radio power is 2 
orders of magnitude greater than that of typical Seyferts. 
The mid IR image is dominated by an unresolved core. 
The ISOPHT spectrum shows weak PAH emission. The 
TIMMI2 spectrum is badly affected by telluric sam- 
pling noise and no PAH features are visible. [Ne II] and 



[Ar III]/[Mg VII] is detected by ISOSWS but not with 
TIMMI2. 

5.2.12. IC4329A 

This Syl galaxy is point-like in our TIMMI2 images, but 
shows extended [O HI] and Ha emission (Mulchaey et 
al. 1996). The TIMMI2 and ISOPHT spectra appear to 
trace the same components, except for the [S IV] line 
which is stronger in the large aperture (Fig. 115(1 . [Ne II] 
is seen with TIMMI2; the feature near 9.4/im is probably 
telluric. 

5.2.13. Mrk509 

The emission of this Syl galaxy is at 11.9/xm dominated by 
an unresolved core and marginal resolved emission below 
the 10% level of the peak flux (Table 4) is found toward 
the South-East (Fig. This may be supported by the 
slight increase of the mid-infrared emission in the larger 
ISOPHT aperture as compared to the TIMMI2 spectrum 
(Fig. El- Source extension towards the south is also re- 
ported by near infrared observations (Winge et al. 2000). 
The ISOPHT spectrum is typical for a Syl (Clavel et al. 
2000), it shows weak PAH emission over a strong contin- 
uum. In the TIMMI2 spectrum PAH bands are not de- 
tected (Table 8). Both spectra do not show evidence of 
silicate absorption. [Ne II] is detected by ISOSWS but 
not with TIMMI2. 

5.2.14. NGC6240 

A famous ultraluminous merger with a deeply hidden 
AGN. The AGN reveals itself by the strong iron line at 
6.4 keV and the powerful continuum at 100 keV (Matt et 
al. 2000). Two hard X-ray nuclei have been recently de- 
tected by Komossa et al. (2003). Near infrared spectra 
indicate very low metalicity {^1/10 solar, Oliva & Origlia 
1998). NICMOS images (Scoville et al. 2000) show two nu- 
clei separated by 1.6" (0.8 kpc) North-South. There are 
also two nuclei at radio wavelengths 1.4" apart (Condon 
et al. 1996), the CO (2-1) peak lies in between (Tacconi 
et al. 1999). The TIMMI2 image at 11.9 /im (Fig. 8) re- 
veals the nuclear region to be dominated by a strong core, 
which we identify with the Southern nucleus (Tecza et al. 
2000). A fainter structure, at 10% of the peak flux of the 
core, extends ^2" to the North and matches the structure 
observed at lower resolution at near-infrared wavelengths 
(Tecza et al. 2000). 

The ISOPHT spectrum is dominated by strong PAH 
emission and shows in addition emission lines of [Ar II], 
Ha 0-0 S(5) (blended with [Ar II]) and H2 0-0 S(3). 
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Table 8. Continuum subtracted 11.3/im PAH band fluxes, 
F, (in 10~^°W/cm^) and PAH band-to-continuum ratios, 
i?, of the TIMMI2 (3" slit) and ISOPHT (24" aperture) 
spectra. 
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Interestingly, the latter line is absent in our TIMMI2 spec- 
trum. Most of the [Ne II] emission (~2/3) originates from 
outside the 3" slit (see Table 7). Both spectra (Fig. [TK|l 
show prominent PAH emission. The 8.6/im/7.7^m band 
ratio is small in ISOPHT and the PAH band at 8.6/im 
feature is not even detected with TIMMI2. Between 9 and 
10.5/xm the TIMMI2 spectrum is noise dominated. The 
sihcate absorption, as of ISOPHT suggests Ay > 25 mag. 

5.2.15. NGC7674 

This Sy2 is slightly extended in the radio and has within 
0.5" a double source (Kukula et al. 1995). The decon- 
volved 11.9/im TIMMI2 image is also extended and at 0.5" 
(Fig. Table 5) dominated by a single core. Previous 
mid IR maps by Miles et al. (1996) have insufficient res- 
olution to resolve the galaxy but their absolute photome- 
try is consistent with ours. The 5" resolution UCL spec- 
trum (Roche et al. 1991) is rather featureless, whereas the 
ISOPHT spectrum contains several PAH bands (Schulz, 
priv. communication). We could not yet perform TIMMI2 
spectroscopy. 

6. Heating and evaporation of PAHs 

We discuss the survival of PAHs in starburst and AGN 
environments with energetic photons of, at least, a few 



Fig. 16. The absorption efficiency of a graphite sphere of 
lOA radius. Mie theory is still assumed to be valid at X- 
rays. Optical constants from Laor & Draine (1993). 

hundred eV. Such photons can only be emitted by an X- 
ray source. 

6.1. The evaporation condition 

Let P{T) be the temperature distribution of an ensemble 
of very small grains, for example, PAHs, so that P{T) dT 
gives the probability of finding a particular grain in the 
temperature interval dT around T. Whereas grain cooling 
usually proceeds through emission of infrared photons, in 
a very hot particle, above '--^2500K, evaporation is more 
effective (see Fig. 8 of Voit, 1991). The evaporation rate 
is proportional to the number of surface atoms, A^surf , for 
which approximately A^surf ~ N^^^ if N is the total num- 
ber of atoms in the grain. 

Let Pi and P2 be the probabilities that a surface 
atom is at an energy level below or above the binding 
energy E^,, respectively. Usually, P2 <^C 1, so Pi ~ 1, 
and obviously Pi + P2 = 1- For a Boltzmann distribu- 
tion, P2/P1 ~ exp{— / kT) . Due to phonon collisions, 
an atom changes its level at a rate given by the Debye fre- 
quency t'D- In one second, the atom is P2i^B times above 
-Bb, the evaporation rate is therefore 

Rev = Nsnrf vj, e-^^l'^'^ . (1) 

Evaporation rises exponentially with temperature. To 
counterbalance it, atoms must be built into the grain. The 
rate at which this occurs is proportional to the thermal 
velocity of atoms, or to T^^^, and to the density of the in- 
terstellar medium. The exponential process of evaporation 
wins over accretion by orders of magnitude in a fairly nar- 
row range of T, the density is not relevant. Let us assume 
that the grains evaporate if a certain fraction of them, /ev, 
is above some critical temperature Tcr- The latter should 
be somewhat greater than Tov . Relegating all complicated 
or unknown physics to the parameters /ev and Tcr, we 
write the condition for evaporation as 

/•OO 

/cv = / P{T) dT > /ev . (2) 

Reasonable numbers for Ta in Q are 2500 K and /cv ^ 
10^^, but such suggestions do, of course, not constitute a 
theory. 

6.2. X-ray cross sections for dust absorption 

An AGN differs from a starburst nucleus of the same bolo- 
metric luminosity by the hardness of its radiation field. 
We are here interested in the question how the presence 
of very energetic photons in AGNs affects the chance of 
survival of small grains. First, we remark briefly on the 
X-ray absorption cross section of dust. 
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Fig. 17. The distribution function of the temperature, 
P(T), for graphite grains of 5A radius. Monochromatic 
flux F — 10** erg cm~^ s^^, photon energies from 10 to 
80eV. 



When a grain absorbs an optical or UV photon of fre- 
quency v, its internal energy rises by the full amount hv. 
Furthermore, when the grain is not smaller than the wave- 
length, the absorption efficiency is approximately one. The 
situation is different at X-rays . The interaction of hard 
photons with interstellar dust and the limits of Mie theory 
have been studied by Voit (1991), Laor & Draine (1993) 
and Dwek & Smith (1996). Below lOkeV, the main ab- 
sorption process is photo-ionization of inner atomic shells. 
The excited electron loses energy in inelastic scattering as 
it travels through the grain, but may also leave the grain 
and carry away kinetic energy as a photo-electron. When 
an X-ray photon has created a gap in the innermost shell, 
the gap will be filled by a downward transition of an elec- 
tron from an upper shell and the energy may escape ei- 
ther as a photon, or non-radiatively through ejection of 
an Auger electron. 

The computations of Dwek & Smith (1996) indicate 
that for graphite particles of 50A radius, soft X-ray pho- 
tons with hf < 100 eV deposit practically all their energy 
in the grain and Mie theory is still valid. For photon en- 
ergies above hv = 100 eV and particle radii a < 50A, only 
part of the photon energy is deposited in the grain and 
we therefore reduce the absorption efficiency by a factor 
C!C l/v. 

Fig. 1 161 illustrates the absorption efficiency, Q'^^^, from 
Mie theory of graphite grains with radii from 10 to lOOOA 
(see Fig. 2 of Laor & Draine (1993) for a wider range 
of parameters). The jump at 43A (290 eV) is due to the 
ionization threshold of the K-shell of the carbon atoms. 
At still higher frequencies, one has the typical i^'^ cross 
section for ionizing solitary atoms. The drop in Q^^^ at 
wavelengths short of the first ionization level at ^10 eV 
is less steep (Q'^^'^ oc for a < lOOA); for big grains 
(a > lOOOA), Q^^^ is essentially flat at these wavelengths. 

When, as for X-rays, absorption is due to the excita- 
tion of inner atomic shells, the cross section of a grain, 
K^^^, depends only weakly on the way the atoms are 
held together in the solid. One can therefore roughly set 
^abs proportional to the number of atoms in the grain. 
The dust particle in Fig. El with lOA radius consists of 
Nc = 520 carbon atoms and has Q^^^ c± 0.2 in the inter- 
val lOeV <hiy< 20eV, and Q^'^" oc i/"^ ^ for > 20eV. 
For the cross section per carbon atom, C^^^, this im- 
phes C^'^" ~ 10-"cm2 for lOeV < hv < 20eV, and 
^abs ^ ^Q-n X (20/hvf -^c\rL^ for hv > 20eV. We wih 
use these crude numbers to calculate the cross section of 
PAHs. 



Fig. 18. The fraction of graphite grains of SA (solid lines) 
and lOA (dashed lines) radius above a critical tempera- 
ture Tcr = 2500 K (see Equation EJ . The strength of the 
monochromatic flux F is indicated. 

Fig. 19. The fraction of PAHs with iVc = 100 carbon 
atoms above a critical temperature Ta = 2500 K (see 
Equation 121) . The strength of the monochromatic flux F 
is indicated. 

6.3. The temperature distribution of PAHs in a hard 
radiation field 

To investigate how the distribution function P{T) and, 
in particular, the evaporation condition Q depend on the 
photon energy, we consider a source with a monochromatic 
flux of strength F, all photons have then the same energy 
hv. While keeping the value of the flux fixed, we change 
hv, so if iVph is the number of photons, Npi^hv — F stays 
constant. Low energy photons impinge more frequently, 
but raise the enthalpy H of the grain only by a small 
amount. X-ray absorption, on the other hand, is a rela- 
tively rare event, but accompanied by a big jump in H . 
Temperature distributions P{T) are shown in Fig. ll7l for a 
flux F = 10^ erg cm^^ s^^ corresponding to a luminosity 
L = 2.4 X 10^-^ L© at a distance of 100 pc. The photon 
energies vary between 10 and 80eV. Grains much above 
the evaporation temperature of the solid (~2000K) are, 
of course, unrealistic. 

The distributions P{T) in Fig. El display a kink at 
the temperature where the thermal energy of the grain 
equals hv. Still higher grain (momentary) temperatures 
are achieved when a second photon is absorbed before the 
grain has cooled off from the excitation by the first. 

Fig. El displays the integral lev in equation Q as a 
function of photon energy hv for various fluxes and two 
grain sizes. The important point is here that for the small- 
est graphite grains with 5A radius (65 carbon atoms), 
lev grows very rapidly in the range from 20 to 40eV 
and then declines only gradually. The energy of photons, 
{hv), emitted by the hottest O stars (Toff = 5 x 10"* K) 
is only 16 eV when averaged over the total spectrum, and 
equals 22 eV when averaged over the Lyman continuum 
(A < 912A). These numbers are thus upper limits for the 
radiation field in a starburst. Furthermore, the stellar flux 
of an O star plummets at wavelengths below 228A, or 
above hv — 54.1 eV, when helium becomes doubly ion- 
ized. AGNs, on the other hand, have a spectrum typically 
declining with v^^'"^ and emit copious amounts of extreme 
UV and soft X-ray photons as witnessed by the observa- 
tion of fine structure lines from ions with high (> 100 eV) 
ionization potentials. 

We conclude from Fig. El that small graphite grains 
(5A radius) survive in a starburst as long as fluxes are 
below ^ 10'* erg cm^^, but they will be destroyed near an 
AGN of the same flux F. The fate of PAHs of equivalent 
mass is similar (see Fig.EJ- Bigger grains (see curves for 
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Fig. 20. Optically thin dust emission (full line) between 5 
and 15 /zm in a starburst (top) and an AGN (bottom) at 
various distances to the central heating source. The total 
luminosity is L = IO^^Lq. The emission of large grains, 
representing the continuum, is shown dashed. 

lOA radius in Fig.ll8|l are much more resistive. They need 
higher fluxes to get vaporized and survive in a starburst 
environment. 

7. Radiative transfer models 

Here we present comprehensive model calculations for the 
mid-infrared emission of starburst and AGN radiation en- 
vironments in the optical thin and optical thick regime. 

7.1. Optically thin dust emission 

In a first attempt to understand and to illustrate the in- 
fluence of the hardness of the radiation field, we compute 
the optically thin dust emission at various distances from 
the source. A point-like heating source of total luminos- 
ity L = IO^^Lq has either the spectrum of an OB star 
(T* = 25 000 K) to mimic a starburst, or of an AGN with 

oc 1/"° '' in the range from lOA to ~2^m (the exact 
upper limit is unimportant as the dust has little extinc- 
tion there). The distance to the source, and thus the flux 
which the dust receives, is varied from 4 to 100 pc. 

For simplicity we consider in this section only two dust 
components: Large carbon grains with a n(a) oc a~^'^ size 
distribution and radii between 300 and 2400A, and PAHs 
consisting of 100 C and 28 H atoms and an abundance of 
10% of the large carbon grains. The spectra are shown in 
Fig. 1201 They are in arbitrary units and normalised at the 
maximum. 

One notices that large grains become hotter as the 
radiation field hardens. At a distance of 4 pc to the heating 
source, PAHs invariably evaporate, at 10 pc evaporation 
occurs only near an AGN, and at 100 pc they survive in 
both environments, starburst and AGN. 

7.2. Optically thick dust emission 

Nuclei of luminous infrared galaxies are generally dust en- 
shrouded and not transparent. Radiative transfer calcu- 
lations of optically thick dusty nuclei have been carried 
out in various approximations (Kriigel & Tutokov 1978, 
Rowan-Robinson & Crawford 1989, Pier & Krohk 1993, 
Loar & Draine 1993, Granato & Danese 1994, Kriigel 
& Siebenmorgen 1984, Siebenmorgen et al. 1997, Silva 
et al. 1998, Efstathiou et al. 2000, Ruiz et al. 2001, 
Siebenmorgen et al. 2001, Nenkova et al. 2002). 

In the subsequent spherical models, for starbursts and 
AGNs, the total luminosity is again 1O^^L0. The dust 
cloud has a constant density, a visual extinction of Ay 
25 mag (from the surface to the centre), and an outer cloud 




6 8 10 12 14 16 
wavelength (|j,m) 

Fig. 21. Comparison of radiative transfer model spectra 
of dusty starbursts (solid) and AGN (dashed). In both 
models the visual extinction towards the cloud centre is 
= 25 mag and the luminosity is i = lO^^L©. Spectra 
are calculated for a distance of 50Mpc. 




1 10 100 

radius (pc) 

Fig. 22. Normalised intensities at 12 (solid lines) and 
25/im (dashed lines) for starbursts and AGN models as 
presented in Fig. [21 At both wavelengths the predicted 
FWHM of the starburst galaxy is much larger than for the 
AGN. 



radius of 500 pc. To be a bit more realistic the dust con- 
sists here, besides large carbon grains, also of large silicate 
grains of the same size distribution and with a 12% larger 
total volume. Absorption and scattering efficiencies are 
calculated using Mie theory and optical constants by Loar 
& Draine (1993). In addition, there are small graphites 
and small silicates of 10 A radius, with an abundance of 
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0.5% that of the large carbon grains. PAH parameters are 
as described in the previous section. 

The radiative transfer in a dust cloud around a star- 
burst is computed following Siebenmorgen et al. (2001), 
around an AGN after Siebenmorgen et al. (1997) and the 
inclusion of soft X-rays as in this work. The major differ- 
ence between both radiative transfer models is that AGN 
are heated by a central source whereas starburst galaxies 
are heated by stars which are distributed over a large vol- 
ume. In case of a starburst, each of the OB stars have a 
uniform luminosity of 20 000 L0 and a stellar temperature 
T* = 25 000 K. The stellar density changes with galactic 
radius r up to 500 pc from the center like oc 1/r (there 
must be a cutoff at r = 0, but its exact position is unim- 
portant). Each OB star is surrounded by a hot spot of 
constant dust density, corresponding to n(H) 10*cni~^. 
The size of the hot spot is determined by the condition 
of equal heating of the dust from the star and from the 
radiation field in the galactic nucleus. As for the AGN 
models, the inner radius of the hot spots is given by the 
photo-destruction or evaporation radius of the grains. 

Our model spectrum for the starburst and the AGN 
is displayed in Fig. |^ The starburst shows strong PAH 
emission and a continuum which rises steeply beyond 
10/im, whereas the AGN is almost free of PAH emission. 
Cross cuts over the model nuclei, at various wavelengths, 
are shown in Fig. At 12/im and 25/im, the AGN is 
much more compact than the starburst. Therefore the 
mid-infrared size may be used to distinguish between the 
two kinds of energy sources. This statement holds unless 
the stars are distributed by more than a few pc. 

8. Conclusion 

We presented mid infrared imaging and spectroscopy for 
a sample of 23 galactic nuclei. The observations were 
obtained with the ISO (14" - 24" resolution) and with 
TIMMI2 at the 3.6m ESO telescope (3" and 0.5" reso- 
lution for spectroscopy and imaging, respectively). The 
galaxies are mostly starbursts and AGNs, with apparent 
mid-infrared fluxes of 20 mJy or more. Our main results 
are: 

a) Seyfert nuclei are dominated by a compact core 
while starbursts do not display such a central concentra- 
tion; they are usually spatially resolved by TIMMI2. 

b) Almost all our ISO spectra (18 out of 19) are marked 
by bands at 6.2, 7.7, 8.6 and 11.3^m, generally attributed 
to PAHs. In the TIMMI2 spectra, which refer to the in- 
nermost nuclear regions, these bands are present in star- 
bursts, but not detected or significantly reduced in AGNs. 
This corroborates previous suggestions that PAHs evap- 
orate much easier in the hard photon environment of an 
AGN, but are often able to survive in a starburst. 

c) Applying a simple criterion for dust evaporation, 
we quantify the conditions under which small grains are 
photo-destructed in a monoenergetic hard radiation field. 
Analyzing their temperature distributions, wc find that 
PAHs or very small graphite particles with less than 100 



carbon atoms cannot survive if the photons have energies 
above 20 eV and their flux is greater than 10^ erg cm~^ 
s~^. The detection or non-detection of PAH bands can 
thus serve as a diagnosis whether a galaxy is predomi- 
nantly powered by an AGN or a starburst. Of course a 
very dense stellar cluster provides a high radiative density 
and destroys PAH (Siebenmorgen 1993). Therefore the 
dividing line between AGN and starbursts is not sharp, 
as demonstrated by the nuclear spectrum of M83 and 
Circinus. 

d) In radiative transfer models that include soft X- 
rays and incorporate the above dust evaporation criterion, 
the starburst galaxies show PAH bands everywhere, and 
in the mid-infrared they are stronger and more extended 
than AGNs of the same luminosity. 

e) The [S IV] and [Ne II] fine structure lines at 10.5 
and 12.8/im, which can still be excited in an O star envi- 
ronment, are in both types of galaxies generally stronger 
in the large (ISO) aperture spectra indicating that their 
emission comes from outside the nucleus. 
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